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Abstract 

We present a general approach to incorporate hadronic as well as quark degrees of freedom in 
a unified approach. This approach implements the correct degrees of freedom at high as well as 
low temperatures and densities. An effective Polyakov loop field serves as the order parameter 
for deconfinement. We employ a well-tested hadronic flavor-SU(3) model based on a chirally 
symmetric formulation that reproduces properties of ground state nuclear matter and yields 
good descriptions of nuclei and hypernuclei. Excluded volume effects simulating the finite size 
of the hadrons drive the transition to quarks at high temperatures and densities. We study the 
phase structure of the model and the transition to the quark gluon plasma and compare results 
to lattice gauge calculations. 
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1. Introduction 

A central task of ultra-relativistic heavy-ion collisions is to investigate phase transi- 
tions of the strongly interacting matter that is created in the fireball of the collision. 
Experimentally the investigated parameters of temperature and chemical potential at 
which possible phase transitions are probed, can be varied by studying heavy-ion colli- 
sions at different beam energies, from LHC and RHIC energies with very low net baryon 
density to lower energies that will be especially investigated at the FAIR facility at GSI 
which will probe higher densities or chemical potentials, respectively. The most important 
phase transitions occurring in hot and dense matter are the restoration of chiral symme- 
try and the deconfinement transition. In order to model these transitions theoretically 
the main problem arises from the very different degrees of freedom in the limit of low and 
high temperature/density. At high excitation energy the system is described in terms of 
quarks and gluons, whereas at low excitation (or as limiting cases, in the vacuum or in 
the nuclear matter ground state) the effective degrees of freedom are hadrons. In order 
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to describe non-trivial phase transition behavior it is therefore necessary to develop a 
model that contains both sets of degrees of freedom with the correct asymptotics. In this 
article we present a model of that kind, the Hadron-Quark- Model (HQM). We investigate 
the phase transition and present comparisons to lattice gauge calculations at vanishing 
chemical potential. 



2. The HQM model 



The underlying hadronic SU(3) model has the following structure (see [1] for details). 
In mean field approximation one has 

L ^kin ~l~ ^int ~t~ ^rneson ■ (1) 

L int contains the interactions of the baryons and meson fields: 

Lm = ^[70(5^ + + m*)ipi. (2) 

i 

The effective baryon masses m* read 

m* = g ia <7 + g K C + Srrii (3) 

including couplings to the scalar field plus a small explicit mass term. L mesons includes the 
mesonic self-interactions, which in the case of the scalar fields generate the spontaneous 
chiral symmetry breaking, as well as self-interactions of vector fields and an explicitly 
chiral symmetry breaking term: 

T 1 , 2 2, 2 ,2n ( 4 , 4> 4 . o 2 ,2 . ^4> . 2w0 3 



+ \k (^ + C 2 ) - fcx(- 2 + C 2 ) 2 - *2 ( y + C 4 ) - h^C 
+ mlfrra + (^/2m 2 k f k - -^=to 2 ( 



X 4 X 4 cr 2 C 

""J — ^4 j In 2 

Xo Xo ^oCo 



+ X 4 -Xo + ln^4-fc 4 ^ln^f . (4) 



The fields a and ( denote to the non-strange and strange scalar quark condensates, 
and u>, (j> are the corresponding vector fields. The dilaton field x represents the gluon 
condensate in the system. 

In addition the model contains quark degrees of freedom that couple linearly to the 
mean fields together with a Polyakov loop $ field that serves as the order parameter for 
deconfinement in the spirit of the PNJL model for quarks. We adopt a standard choice 
of an effective potential for the Polyakov loop [8] : 

U = -ia(T)$$* + b(T)ln[l - 6$$* + 4($ 3 $* 3 ) - 3($$*) 2 ] (5) 

with a(T) = a T 4 + aiT T 3 + a 2 T^T 2 , b(T) = & 3 T 3 T where the constants arc fitted to 
reproduce quenched lattice results. 
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Fig. 1. Total particle number densities for the different particle species divided by T 3 as a function of 
T at fig = [2]. The black line shows the total number of quarks+antiquarks per volume while the 
green (dotted) line refers to the total meson density and the red (dashed) line to the number density of 
hadronic baryons+antibaryons. 

Thus effectively we couple a hadronic model with a PNJL model for quarks . In order 
to naturally suppress hadrons at high densities and temperatures we take into account 
excluded volume effects in a thermodynamically consistent manner as described in [2,5]. 
For an alternative way to formulate the HQM model and to suppress hadrons at high 
temperatures see [6]. 

3. Results 

We solve the model equations by minimizing the grand canonical potential for given 
temperature and baryochemical potential. Figure 1 shows the resulting value of the par- 
ticle plus antiparticle densities as function of temperature for the case of vanishing chem- 
ical potential. The critical temperature, defined as the maximum of the derivative of the 
scalar field, has a value of T c = 183MeV. One can see that quarks start to dominate the 
system at around T c . However, a quite broad range of temperatures can be observed, 
where the matter consists of a mixture of quarks and hadronic degrees of freedom. This 
is a rather natural outcome given the smooth cross-over transition resulting from this 
model calculation as well as lattice simulations [3,4]. Thermodynamical quantities like 
the interaction measure also compare well to lattice results [2] . 

Studying the phase diagram as function not only of T but also of the chemical potential 
yields results shown in Fig. 2. For the parameters used in this study the corresponding 
nuclear ground state does not contain quarks but is purely governed by the hadronic 
chiral Lagrangian. The phase transition is a cross over for all values (this might change, 
however, if one includes more hadronic resonances). One can see the transition to chiral 
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Fig. 2. Contour plot of the chiral condensate normalized to the vacuum value as a function of temperature 
and quark chemical potential including quark repulsive interactions using a value of the interaction 
strength of g quJ =3 [2]. The dashed grey line indicates where the change of the chiral condensate with 
respect to T and fi q has a maximum. The dashed black line shows the maximum of change of the 
Polyakov loop. Here, there is a well-defined nuclear ground state without quarks. The corresponding 
liquid gas phase transition is shown as the black solid line with a critical endpoint at T cep r& 16 MeV. 

restoration as well as a first order liquid-gas phase transition which continues after a 
temperature of T — 16 MeV as a cross over and which joins the other chiral transition at 
lower chemical potential and higher temperatures. Also shown in the figure is the phase 
transition line to deconfinement, which stays at higher temperatures for large chemical 
potential. 

Various studies of simulations of heavy-ion collisions implementing the equation of 
state of the HQM model have been performed [2,9] that show the importance of the 
quark phase, for instance in the case of dilepton production. More calculations in this 
direction are in progress. 
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